Traces of Thermalization from pt Fluctuations in Nuclear Collisions 



S. Gavin 

Physics and Astronomy Department, Wayne State University, 666 W Hancock, Detroit, MI 48201 

(Dated: February 8, 2008) 

Scattering of particles produced in high energy nuclear collisions can wrestle the system into a 
state near local thermal equilibrium. I illustrate how measurements of the centrality dependence of 
the mean transverse momentum and its fluctuations can exhibit this thermalization. 
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Fluctuations of the net transverse momentum have re- 
cently been measured, with the STAR, PHENIX, NA49, 
and CERES experiments reporting substantial dynamic 
contributions |j, yj 14 ■ Such fluctuations can provide 
information on collision dynamics and, perhaps, the QCD 
phase transition Preliminary PHENIX and STAR 

data in Au-I-Au collisions show that pt fluctuations in- 
crease as centrality increases 0, Importantly, data 
from these same experiments exhibit a strikingly simi- 
lar increase in the mean transverse momentum (pt), a 
quantity unaffected by fluctuations 7, 8]. 

I ask whether the approach to local thermal equilib- 
rium can explain the similar centrality dependence of (pt) 
and Pt fluctuations. My focus is on fluctuations, to de- 
velop the appropriate theoretical tools and experimental 
observables. Dynamic fluctuations are characterized by 
the observable {SptiSpt2) analyzed by STAR ^J], where it 

is termed '^l^p^^ dynam^ 

and CERES For particles of 
momenta pi and p2, one defines 



{5ptiSpt2) = / dpidp2 



P2(Pl,P2) 

{N{N~1)) 



SptiSp 



t2, 



(1) 



where Spti = Pu — (Pt), (■ ■ ■) is the average over events, 
and dp = dydPpt- This definition exploits the relation of 
event-by-event fluctuations to inclusive correlation func- 
tions discussed in The pair distribution is 



P2(Pi,P2) = diV/dpidp2, 



(2) 



where //92c;Pic!p2 = {N{N - 1)) for multiplicity N Q. 
Observe that ^ depends only on the two-body correla- 
tion function 



''(Pl,P2) = /02(Pl,P2) - /3l(Pl)/Ol(P2) 



(3) 



with /3i(p) = dN/dp, since the integral over pipi van- 
ishes due to the definition of Spt- Alternative fluctuation 
observables $pj, Fp^, and Aap^ proposed in 0], 0, and 
|3| measure many-body correlations of all orders. These 
quantities are roughly equivalent 



<^>pja^ A<jpJ<j^ N{Spti6pt2)/2a^ (4) 



when dynamic fluctuations are small compared to statis- 
tical fluctuations a'^ ^ {pD - (pt)"^ |1E3- 



STAR measurements of pt fluctuations in fig. ^ show 
an increase for low multiplicities corresponding to pe- 
ripheral coUisions at s^/^ = 130 GeV 1|. PHENIX mea- 
surements at 200 GeV in fig. also show such an in- 
crease for Fp^ The increase appears to peak and 
possibly saturate for multiplicities corresponding to mid- 
peripheral impact parameters. In addition, the data may 
show a decrease for Act^j and Fp^ for the most central col- 
lisions. While these measurements are preliminary and 
bear large uncertainties, this centrality dependence has 
already been attributed to phenomena associated with 
the QCD transition [Tlllli|. 

I attribute the trend in fig. ^ to the onset of thermal- 
ization in increasingly central collisions, motivated by a 
similar behavior of the measured {pt) in fig.|21,'dl3- Ther- 
malization occurs as scattering between particles pro- 
duced in the collision drives the system toward local ther- 
mal equilibrium. The system is characterized by a phase 
space density /(x, p, t) that varies from collision event to 
event. As the system approaches local equilibrium the 
event-averaged (/) tends toward the Boltzmann-like dis- 
tribution {f^) that varies in spacetime through the tem- 
perature T(x, t). I show here that thermalization alters 
the average transverse momentum following 



{pt)^{pt),S+{pt)e{l-S), 



(5) 



where S is the probability that a particle escapes the col- 
lision volume without scattering. Dynamic fiuctuations 
depend on two-body correlations and, correspondingly, 
are described by 

{6ptiSpt2) = {5pti5pt2) oS"^ + ((5pfi'5pt2>e(l - S)'^. (6) 

The initial quantities {pt)o and {5ptiSpt2)o are deter- 
mined by the particle production mechanism, while {ptje 
and {5pti5pt2) e depend on the state of the system near 
local equilibrium. 

To understand how thermalization can cause the com- 
mon trends in figs. Q and [21 observe that as central- 
ity is increased, the system lifetime increases, eventually 
to a point where local equilibrium is reached. Conse- 
quently, the survival probability S* in Q and ^ de- 
creases from unity as the impact parameter decreases. 
Both (0) and © peak for impact parameters near the 



2 



point where equilibrium is established. The behavior in 
collisions at centralities beyond that point depends on 
how subsequent hydrodynamic evolution changes {pt)e 
and {SptiSpt2)e as the system size and lifetime increase. 
Systems formed in the most central collisions can expe- 
rience cooling that reduces © and ©. 

For both the average pt and its fluctuations to increase 
during thermalization as in figs. ^ and |21 both {pt)e and 
{SptiSpt2)e must exceed the initial values. For the av- 
erage transverse momentum, this implies that the tem- 
perature T at thermalization must be quite high, since 
{pt)e ocT. A value {pt)o ~ 350 MeV near that measured 
in pp collisions implies T ~ 400 MeV, suggesting that 
partons contribute to thermalization. 

In the following paragraphs I estimate {5pti5pt2) o and 
{5pti5pt2)e- Next, I formulate a nonequilibrium approach 
capable of treating fluctuations based on the Boltzmann- 
Langevin equation in the relaxation-time approximation. 
Here, I sketch the derivation of Q and jnj, leaving the 
details for a longer paper. 

Transverse momentum and particle density fluctua- 
tions arise partly due to the particle production mecha- 
nism, e.g., string fragmentation. These fluctuations were 
measured in proton-proton (pp) collisions 0. To use 
these pp results to estimate {5pti5pt2)o for nuclear col- 
lisions, I apply the wounded nucleon model to describe 
the soft production that dominates {pt) and (SptiSpt2)- 
The charged particle multiplicity N and other extensive 
quantities are assumed to scale linearly with the number 
of participant nucleons M, while the intensive one-body 
observable (pt) is independent of M. Centrality is deter- 
mined by iV/iVmax ~ M{b)/M{0) for impact parameter 
b, averaged over collision geometry. 

To estimate the initial {SptiSpt2) using the wounded 
nucleon model, I follow the appendix in ref. |^ to obtain 



{SptlSpt2)o = 



2{SptiSpt2] 



pp 



M 



1 + R 



■pp 



1 + Raa J 



(7) 



The term outside the parentheses is expected because 
(|H) measures relative fluctuations and, therefore, should 
scale as M~^; note that pp collisions have two partici- 
pants. The term in parentheses accounts for the normal- 
ization of JU to (N(N - 1)) = (iV)2(l + Raa) rather 
than (N)^. From |9|, the robust variance Raa satisfies 



Raa 



dpidp; 



r(pi,p2) {N^)~{Ny-{N) 



(8) 



and scales as Raa oc M 
ply {SptiSpt2)pp/{p, 



'PP 



ISR measurements im- 
0.015 [13. HIJING gives 
Rpp ^ 0.45 and Raa ^ 0.0037 for central Au-|-Au for the 
rapidity interval Arj — 1.5 studied in 0. To compare 
Q to N (6ptiSpt2) / {pt)^ in fig. I assume central col- 
lisions produce N « 825 charged particles in Arj = 1.5, 
i.e., dN/dr] « 550. 

Near local thermal equilibrium, dynamic fluctuations 
occur because initial state fluctuations result in transient 



spatial inhomogeneity that can survive thermalization. 
The inhomogeneity would eventually disappear due to 
diffusion and viscosity, but can be observed if freeze out is 
sufficiently rapid. Inhomogeneity is essential for dynamic 
fluctuations, since {SptiSpt2) and ^p^ would otherwise 
vanish for p2 = PiPi- 

To see how inhomogeneity can survive thermalization, 
observe that local equilibrium is achieved when the av- 
erage phase space distribution of particles within a small 
fluid cell (/) relaxes to the local equilibrium form (f^). 
The time scale for this process is the relaxation time 
discussed later. In contrast, density differences between 
cells must be dispersed by transport from cell to cell. The 
time needed for diffusion to disperse a dense ffuid mass 
of size L ^ (|Vn|/n)^-'^ is ^ i/L^/u^j^, where v^^ ~ 1 
is the thermal speed of particles. This time can be much 
larger than z/~^ for a sufficiently large fluid mass. The 
rapid expansion of the collision system further prevents 
inhomogeneity from being dispersed prior to freeze out. 

Inhomogeneity produces spatial correlations: it is more 
likely to find particles together near a dense fluid mass. 
These spatial correlations entirely determine the phase- 
space correlations when the momentum distribution at 
each point is thermal. I write 

r{Pi,P2)^ y"dxidx27'(xi,pi,X2,P2,t) (9) 

evaluated at the freeze out proper time Tp, where the 
phase-space correlation function is 



^12=(/l/2)-(/l)(/2)-<5l2(/l), 



(10) 



for (5i2 = (5(xi — X2)(5(pi — P2). A small change in density 
Sn will initially drive the system from equilibrium by an 
amount Sf^ — f^Sn/n. The corresponding phase space 
correlations are described near equilibrium by 



^f2 = MM'^(xi,X2), 



(11) 



(ni) (712) 

where the spatial correlation function is 

r(xi,X2) = (^1^2) - {ni){n2) - S^ini). (12) 

The form (|ll|l ensures that both 'Pf2 and r(xi, X2) vanish 
in global equilibrium, where particle number fluctuations 
obey Poisson statistics. I use and (^1- 1)12(1 to find 



{SptiSpt2)e= y"dxidx2 r(xi,X2) 



(5pt(xi)(?pt(x2) 

{N{N-1)) 



(13) 



where the local transverse momentum excess, Spt{x) = 
J dp {pt — (pt))/(x, p)/n(x), vanishes if the collision vol- 
ume is uniform. 

To estimate {6ptiSpt2)e using (|13|l . I assume that 
Bjorken scaling holds and that longitudinal and trans- 
verse degrees of freedom are independent. I then write 
the transverse coordinate dependence of p2(l as 

r(xi,X2) cx girti)g{rt2)ci\rti - rtzl) (14) 
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near equilibrium 
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FIG. 1: (a) Dynamic pt fluctuations computed using com- 
pared to STAR data Q- (b) Same for PHENIX data ^. 



where the density is n(xi) cx g(rt). I parameterize g and 
c to be Gaussian with r.m.s. widths Rt and ^, respec- 
tively the transverse radius and correlation length. The 
momentum excess Spt{'x) in depends on the temper- 
ature profile of the system, since J ptf{x,p)dp/n{x.) oc 
T{rt). Similarly, (pt) (x ||r||, for the density-weighted 
average ||r|| = / g{rt)T{rt)drt, so that 



Sptin) = {pt)[f{rt) 



(15) 



I parameterize T{rt) = r(rt)/||T|| as Gaussian of width 
Rp and use n oc to fix Rp = VSRt. 

The dynamic pt fluctuations near local equilibrium 
then satisfy 



{Spti5pt2) 



F 



{pt?RAA 

1 + Raa 



(16) 



where Raa is given by ((SJ . The quantity F is dimension- 
less and depends on the ratio of the correlation length 
to the transverse size Rt ■ I use (|14f) and (|15|l to compute 
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FIG. 2: Average pt from (|5J| compared to data 



approximating the collision term /[ f ] using a single re- 
laxation time v^"^ . Following 0,^3i I longitudinal 
boost invariance to write the left side of (|18|l as df j dr at 
fixed PzT. Longitudinal expansion further implies that 
the density satisfies n (r) oc t~^, while {pt)e oc T oc t~'^ 
for < 7 < 1/3; see [13 • I then multiply both sides of 
(18) by |p(| and integrate over momentum to obtain 



{Pt) = {Pt)oS + 

ct - 

The survival probability is 



<Pt)l 



S 



-r 



{ro/TFr, 



(19) 



(20) 



where u = {av,-c\)n{T), a = i'oTq for the formation time 
To, the scattering cross section is a, and v^d is the relative 
velocity. For relevant values a 3> 7, I approximate H19() 



by with (pt 



{ptmro/rr- 



To compute the evolution of {SptiSptT), I obtain re- 
laxation equations for V . Fluctuations due to scattering 
and drift are described by adding a Langevin force to the 
right side of (|18|l 15]. On a discrete phase space lattice 
Pi , Xj , the Boltzmann-Langevin equation is 



dh/dT = -uu, - in + Q 



(21) 



F = 



where Ci(''") is a Langevin force. To incorporate the effect 
of fluctuations near local equilibrium, I further treat /f 
as a stochastic variable subject to an additional Langevin 
force, so that dff/dr = Xi, plus a diffusive relaxation 
term that I need not specify for a diffusion time scale 
td ^ i^^^. The Langevin terms satisfy {Ci{''')Cj {'t')) — 



i\rti - rt2|)[T(rti) - l][T(rt2) - 1]||, (17) + r)5y(5(T - r') and (G(t)Xj(t')) = -i^rS,jd{ 



a double density- weighted average over rti and rt2- I 
find F = 0.046 for ^t/Rt = 1/6. To determine ^ for 
fig. Uli, I take Raa = 0.0037 and N « 825 as before. I 
emphasize that the HIJING Raa value builds in fluc- 
tuations from resonance decay and, moreover, is roughly 
consistent with measured net charge fluctuations ^14]. 

Let us now describe the fast local relaxation of the 
phase space density f to f^. I start with a Boltzmann- 
like kinetic equation 



r'), as required by detailed balance for the relaxation- 
time collision term .15j . The Boltzmann equation used 
to compute (pt) is the mean value of (|21|l . 

I use standard methods to obtain the following 
two-body relaxation equations 



dVij/dr 
dCij/dr 



(22) 
(23) 



a//9i + vp.v/ = /[/]«-i.(/-r). 



(18) 



where I introduce the auxiliary function Cij = (fifj) — 
{fi){fj)- Observe that Vij = Cij = in global equi- 
librium where the time derivatives vanish. I solve (22) 
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and (23) assuming that Cy initially vanishes and obtain 
{SptiSpt2) from (1), (3), and (9). Equation 10) follows, 
but is exact only if one neglects the time dependence of 
(|16|l implied by {pt)e oc t~'' . For a ^ 7, I approximate 
this dependence in © by taking {SptiSpt2)e t"^'*'. 

I now fit this transport framework together with my 
earlier assertion that near-equilibrium pt correlations are 
induced by spatial inhomogeneity, i.e., eq. (11). In the 
relaxation-time approximation is arbitrary, as is f^. 
To deduce either from transport theory, one must use 
(|21|l with the full collision term /[/]. Following yields 
= {f!){fi)0, where 9 = ai2 + E;.. b^i2PiP2. for a and 
6'"^ functions of Xi and X2- In a uniform system these 
coefficients are constant, so that (1), (3) and (9) imply 
{5pti5pt2) = 0, confirming our intuition. Our physically- 
motivated (|ll|l takes 9 « r{xi,X2) /n{xi)n(x2)^ which is 
adequate for our estimate (16). 

Calculations in figs. ^ and [5] illustrate the common ef- 
fect of thermalization on one-body and two-body pt ob- 
servables. Equation © together with the computed (0) 
and (|13|l is in good accord with data. The solid curves in 
all figures are fit to STAR fluctuation data and {pt) data 
(except for N ^ I ignore any energy dependence). I assume 
a = 4 and 7 — 0.15 in central collisions, and parameter- 
ize S{M) by taking a oc M^/^ and Tp - tq (x M^/"^. In 
this work it is not necessary to specify whether the equi- 
librating system is partonic or hadronic. That said, in 
fig- Ell take the same a for all species, as appropriate for 
parton scattering. Measurements of pt fluctuations for 
identified particles can further test whether thermaliza- 
tion is species independent. 

In comparing to PHENIX data in frg.Ut, note that the 
magnitude difference with fig. ^ follows from the differ- 
ent acceptance of STAR and PHENIX. The solid curve 
in fig. agrees with the data within the uncertainty, 
but the dashed curve shows better agreement for 7 — 0.2 
and Tp — Tf) (X M . While agreement with {pt) data for 
the new parameters is less compelling than fig. 12 results 
still fall within the uncertainty. 

Preliminary data from refs. and 0, 111 show tan- 

talizing similarity to the calculations. However, experi- 
mental uncertainty must be reduced to firmly establish 
the low multiplicity rise as well as the behavior at high 
multiplicity. Contributions to {pt) and {SptiSpt2) not in- 
cluded in this exploratory work are diffusion, collective 
radial flow, Bose-Einstein (HBT) correlations, and collec- 
tive hadronization. Collective effects can be important in 
central collisions, where the matter evolves after equili- 
bration. Flow can enhance the fluctuations, while diffu- 
sion can reduce them. HBT effects can be experimentally 
estimated by cutting on each pair's relative momentum. 
This contribution is of order 10% at RHIC energy 
but may be larger at lower energies Q ■ While resonance 



and hard-scattering contributions to fluctuations are es- 
timated by taking Raa in and ((TB|l from HIJING, 
chemical equilibration may modify the centrality depen- 
dence for resonance production, altering {SptiSpt2)- 

Experimental indications that nuclear collisions pro- 
duce matter near local equilibrium are scant and circum- 
stantial. Any experimental evidence of the onset of equi- 
librium — particularly at the parton level — will validate 
those indications. Rapidity dependence measurements 
can distinguish the thermalization effects proposed here 
from alternative explanations [lllll^ . Here, the rapidity 
dependence arises from the dependence of Q and l(T?)| 
on Raa, which is itself measurable 01. 
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